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ABSTRACT 

The  feasibility  of  performing  nondestructive  analysis  of  spent 
power  reactor  fuel  by  means  of  gamma  ray  spectroscopy  was  esta- 
blished.  Dose  rates  to  be  expected  from  spent  power  reactor  fuel 
pins  offered  by  the  AEC  for  study  were  calculated.   Parallel  cal- 
culations were  performed  for  MITR  fuel  elements  and  the  results 
compared  with  measured  dose  rates  of  the  MITR  fuel  elements. 
Gamma  ray  spectra  were  obtained  from  MITR  fuel  elements  and  used 
to  estimate  the  expected  gamma  ray  spectra  for  the  spent  pins. 
Calculations  were  performed  to  determine  the  number  of  Pu-239 
fissions  which  had  occurred  in  the  spent  fuel  pins  and  the  re- 
sults used  to  calculate  the  expected  change  in  the  spectra  due 
to  the  presence  of  Pu-239  fission  product  gamma  rays.   Recommendations 
are  made  to  enhance  the  ability  of  M.I.T.  researchers  to  perform 
the  subject  experiments  on  actual  spent  power  reactor  fuel. 
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I:   INTRODUCTION 

The  possibility  of  utilizing  gamma  ray  spectroscopy  for 
determination  of  fission  product  content  for  use  in  ascertain- 
ing the  irradiation  history  of  spent  fuel  has  been  studied  by 
many  investigators  (1-4) .   Sovka  (5)  utilized  the  excellent 
energy  resolution  of  Ge(Li)  detectors  to  develop  a  method  for 
determining  spatial  distribution  of  absolute  neutron  flux, 
neutron  exposure,  total  U-2  35  burn up  and  irradiation  time  of 
the  fuel.   This  method  was  unique  in  that  it  did  not  require 
absolute  calibration  of  the  detecting  system  since  it  used 
ratios  of  fission  product  activities.   However,  because  it 
does  use  ratios  of  fission,  product  activities,  more  than  one 
fission  product's  gamma  ray(s)  must  be  resolved.   This  work 
was  performed  using  the  spectra  obtained  from  spent  MITR  fuel, 
Since  MITR  fuel  has  a  pre-irradiation  enrichment  of  9  3%,  it 
differs  considerably  from  that  of  the  more  widely  used  and 
economically  more  important  PWR  or  3WR  power  reactor  fuel 
which  usually  has  an  enrichment  of  2-4%.   Because  of  the  im- 
portance of  accurate  knowledge  of  the  present  status  and  past 
history  of  partially  burned  or  spent  power  reactor  fuel  from 
the  viewpoint  of  fuel  cycle  economics,  reactor  core  physics, 
and  most  recently,  safeguards,  it  is  quite  important  that  ad- 
vanced spectrometry  methods  be  extended  to  these  fuel  types. 
This  extension  is  in  part,  one  of  the  aims  of  the  M.I.T.  Re- 
actor Physics  Project  which  is  sponsored  by  the  AEC  and  under 
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whose  aegis  the  present  work  was  performed.   For  example,  the 
power  reactor  fuel  with  its  high  percentage  of  U-238  would 
have  a  significant  number  of  Pu-239  fissions,  and  thus  the 
spectra  from  spent  power  fuel  would  be  altered  because  of  the 
difference  in  the  fission  product  yields  for  U-235  and  Pu-239. 
Self-shielding  effects  upon  the  emitted  spectra  would  also 
obviously  differ.   The  purpose  of  this  study  was  therefore  to 
determine  the  feasibility  of  performing  gamma  ray  spectroscopy 
on  spent  power  reactor  fuel. 

In  order  to  establish  feasibility,  gamma  ray  spectra  were 
obtained  from  spent  MITR  fuel  to  measure  the  energy  resolution 
of  the  detecting  system,  to  obtain  data  from  which  to  estimate 
the  relative  intensities  of  gamma  rays  present  in  power  reactor 
fuel  spectra,  and  to  determine  modifications  to  the  scanning 
system  necessary  for  the  scanning  of  power  reactor  fuel  by 
subsequent  investigators. 

Calculations  of  the  dose  rate  to  be  expected  from  the 
spent  fuel  offered  by  the  AEC  for  study  were  made  in  order  to 
develop  safe  handling  procedures  for  the  highly  radioactive 
fuel,  which  was  in  turn  necessary  in  order  to  obtain  AEC 
licensing  for  the  subject  experiments.   Parallel  calculations 
were  performed  on  spent  MITR  fuel  and  comparisons  made  with 
observed  dose  rates  in  order  to  check,  the  fidelity  of  the  cal- 
culations and  to  normalize  the  predicted  results  to  physically 
measured  quantities.   Finally,  calculations  of  the  relative 

gamma  ray  intensities  present  in  the  power  reactor  fuel  spec- 
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tra  after  a  2.75  year  cooling  period  were  made  and  their  re- 
sults used  to  determine  the  suitability  of  the  proposed  fuel 
for  analysis.   Each  of  the  above  aspects  are  reported  upon  in 
detail  in  the  three  chapters  which  follow. 
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II :   DOSE  RATE  CALCULATION  FOR  SPENT  FUEL 

2. 1  Introduction 

Prior  to  obtaining  spent  power  reactor  fuel  it  was  nec- 
essary to  develop  handling  procedures  and  to  secure  AEC  li- 
censing approval,  which  necessitated  estimation  of  the  dose 
rate  expected  from  the  spent  fuel.   Two  fuel  pins  which  had 
been  tested  in  the  Dresden  power  reactor  were  made  available 
by  the  AEC  for  the  studies.   Some  of  the  known  and  calculated 
information  (6)  on  the  fuel  pins,  numbered  A4  and  A28  in  the 
Dresden  core  loading  is  listed  in  Tables  1  and  2. 

2.2  Background  Information 

Many  studies  (7-11)  have  been  published  on  the  activity 
of  spent  fuel  ranging  from  experimental  studies  based  on  calor- 
imetric  methods  to  theoretical  calculations.   The  experimental 
studies  present  only  total  activities  and  as  such  are  unsatis- 
factory for  dose  rate  claculations ,  since  the  energy  depen- 
dence of  the (activity  must  be  known,  because  the  mass  absorp- 
tion coefficient,  the  buildup  factor  and  flux  required  for  a 
given  dose  rate  are  all  strongly  dependent  on  energy.   The 
theoretical  compilations  have  utilized  fission  product  data 
to  calculate  the  activities  for  varying  irradiation  and  decay 
times.   Blomeke  and  Todd  (12)  have  calculated  the  activities 
of  specific  fission  products,  and  while  this  study  thus  gives 

the  complete  energy  dependence  of  the  spectra,  the  number  of 
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fission  products  makes  this  method  too  cumbersome  for  estimat- 
ing the  dose  rate  from  spent  fuel.   Most  studies  divide  the 
energy  spectrum  of  fission  product  gamma  rays  into  two  or  more 
energy  groups  and  sum  all  gamma  rays  whose  energies  fall  with- 
in these  groups.   If  a  sufficient  number  of  groups  is  used  the 
calculations  will  account  for  the  energy  dependence  mentioned 
previously  without  the  necessity  of  accounting  for  each  fission 
product's  activity  on  an  individual  basis. 

The  study  by  Prawitz  et  al  (7)  was  chosen  to  estimate  the 
activity  of  the  spent  fuel  because  a  sufficient  number  of 
energy  groups  are  used  and  the  decay  of  the  activities  is  pre- 
sented for  cooling  times  up  to  20  years,  whereas  many  studies 
do  not  report  the  decay  for  times  longer  than  one  year.   The 
study  gives  the  activity  of  fission  products  as  a  function  of 
both  irradiation  time  and  cooling  time  and  some  indication  is 
also  given  of  the  trend  for  flux  variations.   This  latter  para- 
meter is  important  only  for  fission  products  which  have  high 
neutron  absorption  cross  sections.   Statistically,  this  group 
of  fission  products  represents  only  a  small  fraction  of  the 
total  fission  product  activity.   If  the  flux  dependence  must 
be  known  the  comprehensive  study  of  Blomeke  and  Todd  (12) 
presents  the  activity  of  spent  fuel  for  various  irradiation 
and  cooling  times  as  a  function  of  neutron  flux.   While  this 
study  does  give  the  dependence  of  activity  on  flux  it  was  felt 
this  advantage  was  offset  by  the  fewer  number  of  energy  groups 
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used. 

2. 3  Method  of  Calculation 

In  order  to  ensure  a  conservative  estimate  and  to  sim- 
plify the  calculations,  it  was  assumed  the  burnup  occurred 
during  a  one  year  irradiation  period,  followed  by  a  cooling 
time  of  700  days. 

The  method  of  calculation  will  be  illustrated  for  fuel 
pin  A28  with  a  4  foot  water  shield. 

1)   The  15  energy  groups  and  their  energy  range  are 
listed  in  columns  1  and  2  of  Table  3.   The  effective  energy  of 
each  group  was  assumed  to  be  equal  to  the  midpoint  energy  of 
each  of  the  groups  and  this  energy  is  listed  in  column  3  of 
Table  3.   This  energy  was  used  for  determining  the  values  of 
the  mass  absorption  coefficient,  buildup  factor  and  the  flux 
to  dose  rate  conversion  factor  to  be  used  with  the  various 
energy  groups.   The  release  of  an  amount  of  energy  equal  to 
the  effective  energy  was  also  assumed  to  represent  one  dis- 
integration.  The  specific  energy  release  rates  for  the  var- 
ious energy  groups  for  an  irradiation  time  of  one  year  and 
700  days  cooling  time  are  listed  in  column  4  of  Table  3.   The 
specific  energy  release  rates  were  converted  to  specific 
activity  by  dividing  by  the  effective  energy  for  each  group 
to  convert  from  Mev/gFP  minute  to  disintegrations/gFP  second. 
This  specific  activity  is  listed  in  column  5  of  Table  3. 
Specific  activity  was  converted  to  Activity  by  multiplying  by 
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the  burnup  experienced  by  the  fuel  pin. 

2)   Once  the  activities  of  the  energy  groups  are  known, 
they  must  be  converted  to  fluxes  present  at  the  point  in  ques- 
tion.  Reference  (13)  contains  expressions  relating  flux  to 
activity,  source  geometry  and  shield  thickness.   For  a  cylin- 
drical source  of  finite  length  the  following  expression  is 
given  for  the  geometry  shown  in   Figure  1: 


<j>  =  B 


Sv^ 


2(a  +  Z) 


F(6,  b2) 


(1) 


where : 

2 

<j>  =  undirected  flux  (y  rays/cm  -sec) 

3 
S   =  volume  source  strength  (y  rays/cm  -sec) 

B  =  buildup  factor  to  correct  for  multiple 

scattering 

r  =  radius  of  cylindrical  source  (cm) 

a  =  distance  from  source  to  observer  (cm) 

z  =  effective  self-absorption  distance  (cm) 

F(<|>,  b0)  =  /e  e"b  sec  6d6 

*     o 

b2  =  bl   +  ^sZ 
b1  =  y1t1 

y .  =  mass  absorption  coefficient  for  shield 

t,  =  thickness  of  shield 

y  -   absorption  coefficient  of  the  source 
Because  the  diameter  of  the  fuel  pins  is  small  the  self-absorp- 
tion distance  is  small  compared  to  the  shield  thickness,  and 

-17- 


o 


o 

CO 

o 

M 


o 

M 

as 


•-J 
o 

03 
E-t 

w 
s 

o 

w 


g 

M 

ft. 


-18- 


the  fuel  pins  may  be  approximated  as  line  sources  of  finite 
length. 

In  this  case,  Equation  (1)  simplifies  to: 


<|>  =  b  -±L 


2-rra 


F(6,   b1) 


(2) 


where : 

ST  =  line  source  strength  (y  rays/cm) 

In  Table  4  the  factors  y1,b1,  B,  F(8;  b  )  ,  BF(9>  b1)  and 

BF(6  b  ) 


are  tabulated  for  a  4  foot  thick  water  shield  for 


2iTa 
each  of  the  energy  groups.   Values  for  \i      and  F(6,b.)  were  ob- 
tained from  Reference  (13) .   Values  from  B  were  obtained  from 
Reference  (14)  for  a  point  isotropic  source  in  water.   The 
value  of  6  used  was  40  . 

In  Table  5  column  2 ,  the  values  of  the  line  source  strengths 
for  each  group  are  tabulated,  having  been  calculated  from  the 
relation: 

c   =  Activity    .  Activity  ,.,> 

L   Pin  Length    40  (2.54) 

The  gamma  fluxes  received  at  the  point  4  feet  from  the 
fuel  pin  were  then  calculated  from  Equation  (2)  and  the  values 
are  contained  in  column  3  of  Table  5.   The  flux  necessary  for 
1  mr/hr  obtained  from  Reference  (15)  is  tabulated  in  column  4 
of  Table  5  for  each  energy  group.   The  dose  in  mr/hr  for  each 

-19- 


Q 

W 
M 

(X 
W 
Eh 


Eh 

Pi 
O 

Pi 
O 
Eh 
U 
<C 

o 

M 

s 

D 

W 
H 
Eh 


W 
PQ 
En 


iH 

ex*  oo 

r> 

r^ 

UD 

U3  UO 

U5 

VD 

VO 

<x>  in 

X! 

1  1 

1 

1 

1 

1   1 

1 

1 

1 

I  I 

^ 

ccj 

w  w 

w 

H 

H 

w  w 

H 

W 

M 

rq  W 

CD 

t= 

r-  rH 

IX) 

ro 

m 

cm  r-- 

r~> 

rH 

in 

r-»  in 

» — 

CN 

in  in 

00 

o 

rH 

H  r- 

ro 

in 

r~- 

^r  rH 

Ua 

m 

rH  r- 

CM 

r^ 

rH 

cn  en 

«* 

m 

ic 

00  rH 

rH 

in 

«* 

•^ 

«3< 

co  cn 

ro 

CO 

CO 

CO  CO 

X! 

u>  I 

1 

1 

1 

l  l 

i 

1 

1 

l  l 

^ 

1  H 

w 

W 

u 

w  w 

w 

m 

rq 

w  w 

CD 

W  <o 

o 

cri 

m 

cn  c^ 

*X> 

co 

CO 

O  CN 

» — 

cn  r> 

CN 

m 

CO 

VD  CO 

CO 

CN 

rH 

in  oo 

Pn 

pq 

H  in 

o 
•  o 

CN 

o 

m 

CO 

rH  CN 

CO 

«nT 

in 

<£>  CO 

o  • 

• 

o 

o 

o  o 

o 

o  o 

o 

o 

^ 

in  m 

O  CN 

rH 

CTl 

VQ 

<y*  r-* 

CN 

• 

• 

•  • 

CQ 

rH  co 

rH 

•^ 

CN 

rH  rH 

rH 

cri 

r-> 

VD  V£> 

•^ 

•^ 

">r 

rH 

CTi 

^D 

m 

LO 

in  i 

1 

i 

•^r 

CO  CO 

XX 

1  m 

1 

i 

! 

i  H 

rq 

w 

1 

1   1 

•» 

W  CN 

rq 

w 

W 

W  o 

o 

o 

rq 

w  w 

CD 

CN   » 

o 

rH 

■^ 

u3  r* 

CO 

r^ 

o 

o  «* 



•  "tf 

pM 

rH  rH 

in 

CN 
O 

rH 

m 

CO 

00  rH 

CM 

^r 

Is- 

rH  rH 

CN  CN 

c^ 

CN 

^r 

rH  ^ 

in 

CO 

o 

in  r- 

*   * 

• 

■ 

rH 

cn  in 

o 

LO 

CN 

oo  m 

rH 

<T>  *# 

rH 

o 

£> 

H  rH 

rH 

rH 

CT\ 

oo  r- 

r- 

KO 

\D 

m  in 

O  O 

LO 

O 

o 

*3<  00 

co 

O 

CO 

O  CO 

00  00 

r- 

"Si" 

m 

f-  h 

Is* 

^r 

o 

oo  in 

in  rH 

CTi 

CO 

r* 

KD   <D 

m 

m 

m 

^  "^ 

rH 

rH  rH 

O 

o 

o 

o  o 

o 

o 

o 

o  o 

P. 

in 

o  o 

o 

o 

o 

o  o 

CN  O 

O 

o 

o 

H  co 

in 

r~ 

CT\ 

rH  CO 

>- 

rH  CO 

in 

r- 

cr\ 

W 

H  rH 

rH 

rH 

rH 

CN  CN 

X 

O  Pi 

(X   D 

W  O  • 

O 

rH  CN 

2«  0 

rH  CN 

n 

«* 

in 

kd  r* 

CO 

c\ 

rH 

rH  rH 

m  o  rs 

-20- 


LO 

w 

PQ 
< 
E-« 


U 

xi 

u 
6 

■*— V 

M 

o 

IT) 

x: 

CM                    k£> 

O 

w  w  \ 

i/iH(NiniriLncNi^Hinu)i^ 

l£> 

W  B   M 

oom     •LO<T>r^ir)'*i-H     »ro 

* 

o  rf  e 

ooiohritno     •rocvjcri"^ 

m 

Q  K  — ' 

II 

W 
H 

H 

g 

x: 

\ 

W 

n 

w 

6 

o 

Q 

H 

romrOCNIfN<N(N(N(N(NCNCN 

A 

II 

WWWWWWWWWWWlxi 

< 

-e- 

OOOOCNOOCOUS^CNOr^ 

Eh 

inHHMn^nnnmnN 

EH 

o 

•                 o            o 

(N        o  o        o           •  o          •  o 

H  oo     •  o  ^     »  ^  v£>     .ooo     • 

•  Or\JU3       •  CN       •  00   <£>       •  O  U3 

>- 

ro     •c^fMt^-oot^-CNi-looc^ro 

■e- 

CNCN»£>r-Hr--rO<Nr-lr-HVOinCN 

o            o 

■H                   rH            00 

Wr~-cy\Wt^wi£icor^t^-oor^ 

(MWW^HOOHWHHHPq 

I^*Xi(N0>iCM00'^'crOCNU3>£i 

-*>^hh^MWHomO 

HNfsiHkDHt^fNOJH^DtN 

x 

O    04 

«  D 

w  o    • 

OHM 

s«o 

HNM^iO^DI^OOOiHHH 

W  O  53 

■21- 


group  is  tabulated  in  column  5  of  Table  5.   The  total  dose  is 
the  sum  of  the  doses  from  each  energy  group. 

Calculations  were  also  performed  for  a  2  foot  thick  water 
shield  and  in  1  meter  of  air.   The  dose  rates  for  fuel  pin  A4 
were  calculated  from  those  of  A28  by  multiplying  by  the  ratio 
of  the  burnup  received.   The  results  of  these  calculations  are 
summarized  in  Table  6. 

In  order  to  check  the  validity  of  the  calculations  a  com- 
parison was  made  between  values  of  dose  rate  calculated  for  2 
MITR  fuel  elements  and  observed  values  for  these  elements. 
The  results  are  tabulated  in  Table  7.   It  was  felt  that  the  dis- 
crepancy was  due  in  part  to  the  approximation  of  the  MITR  ele- 
ment as  a  line  source.   Because  of  the  relatively  large  size 
of  these  elements  appreciable  attenuation  would  occur  in  the 
element.   To  check  this  assumption  the  calculation  was  perform- 
ed again  for  MITR  element  2M15  using  the  finite  cylinder  ex- 
pression, Equation  (1) .   The  radius  was  assumed  equal  to  one 
half  the  thickness  of  the  element,  1.5  inches.   The  calculat- 
ed dose  rate  using  liquation  (1)  was  37.2  mr/hr  compared  to 
43.6  mr/hr  for  the  liir\e  source  expression.   Another  reason 
for  the  discrepancy  between  the  calculated  and  observed 
doses  for  the  MITR  elements  is  the  intermittent  operation 
of  the  M.I.T.  Reactor,  since  Prawitz  (7)  assumed  continuous 
operation  in  the  development  of  the  relations  used  in  these 
calculations . 
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TABLE    6:       CALCULATED    DOSE    RATES    FROM    FUEL    PINS     (MR/HR) 


FUEL 
PIN 

1  METER 
AIR 

2  Ft. 

H2° 

4  Ft. 

H20 

A4 
A28 

12,420 
14,011 

3793 

4280 

38.6 
43.6 
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TABLE  7 :   COMPARISON  OF  CALCULATED  AND  MEASURED  DOSE  RATES 
FOR  MITR  FUEL  ELEMENTS,  4  FT.  HO  SHIELD 


FUEL 
ELEMENT 

TOTALa 
ENERGY 
MWH 

IRRADIATION3 
PERIOD 

COOLING3 
PERIOD 

CALCULATED 
DOSE  RATE 
(mr/hr) 

MEASURED 
DOSE  RATE 
(mr/hr) 

2M15 
2M49 

1296 
1026 

2  yrs. 
1.5  yrs. 

850  days 
130  days 

104 
1455 

26 
450 

(a)  Reference  (2]). 
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Because  of  the  reasonably  good  agreement  between  the 
calculated  and  observed  dose  rate  values  for  the  MITR  fuel  it 
was  felt  the  parallel  calculations  for  the  Dresden  fuel  would 
also  yield  adequate  and  even  more  important,  conservative 
estimates  of  their  expected  dose  rates.   In  particular,  it 
should  be  noted  that  the  expected  dose  rates  from  the  Dresden 
fuel  pins  are  well  below  dose  rates  now  experienced  from  MITR 
fuel  elements  under  similar  shielding  conditions.  This  is  an 
important  conclusion  because  it  implies  that  safe  handling  of 
the  Dresden  fuel  can  be  accomplished  using  available  MITR 
facilities  and  fuel  handling  procedures.   It  is  therefore 
recommended  that  licensing  approval  be  requested  from  the  AEC 
on  this  basis. 

In  view  of  this  conclusion,  it  was  felt  that  no  further 
work  was  required  on  this  aspect  of  the  problem,  and  attention 
could  therefore  be  focussed  on  the  actual  spectrometry,  as  is 
to  be  discussed  in  the  following  chapter. 
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Ill:   EXPERIMENTAL  INVESTIGATIONS 

3.1   Introduction 

In  order  to  use  gamma  ray  spectra  for  analysis  of  spent 
fuel,  certain  fission  product  gamma  rays  must  be  resolvable 
both  in  intensity  and  in  energy  so  that  the  measured  count 
rate  is  a  statistically  significant  measure  of  the  activity  due 
solely  to  the  particular  fission  product.   In  order  to  use  the 
method  of  analysis  developed  by  Sovka  (5)  it  is  desirable  to 
be  able  to  resolve  the  gamma  rays  of  Cs-137,  Cs-134,  Zr-95  and 
Pr-144.   These  fission  products  have  high  fission  yields  and 
relatively  high  energy  gamma  rays;  pertinent  data  concerning 
these  fission  products  is  presented  in  Table  8. 

The  experiments  described  below  were  performed  on  MITR 
spent  fuel  in  order  to: 

1)  Obtain  data  from  which  to  estimate  whether  the  gamma 
rays  for  the  above  fission  products  will  be  resolv- 
able for  cooling  times  on  the  order  of  2  1/2  years 
for  the  power  reactor  spent  fuel  pins  described  in 
Chapter  2. 

2)  Determine  the  improvements  in  resolution  of  the  spent 
fuel  gamma  spectra  resulting  from  the  improved  re- 
solution of  present  Ge(Li)  detectors,  and  use  of  the 
4096  channel  analyzer,  and  thus  ascertain  if  addition- 
al fission  product  gamma  rays  were  suitable  for  analy- 
sis . 
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3)   Check  operation  of  the  scanning  system  and  make  alter- 
ations as  necessary  to  enable  the  equipment  to  scan 
and  hold  spent  power  reactor  fuel  pins. 

3. 2   Experimental  Equipment  and  Procedures 

3.2.1)  MITR  Fuel  Elements 

MITR  fuel  elements  are  of  the  standard  MTR  plate 
type  construction.   A  detailed  description  is  contained  in 
Reference  (  5  ) .   The  important  feature  for  present  purposes 
is  that  the  elements  contain  approximately  160  grams  before 
irradiation  of  highly  enriched  uranium  (93%  U-235) .   Gamma  ray 
spectra  were  recorded  for  two  elements,  7M9  and  5M22.   Some 
of  the  pertinent  information  for  the  elements  are  listed  in 
Table  9. 

3.2.2)  Scanning  Equipment  and  Procedures 

Scanning  equipment  developed  by  Sovka  (5)  for 
scanning  spent  MITR  fuel  was  used  to  obtain  spectra  from 
irradiated  MITR  fuel.   A  schematic  diagram,  shown  in  Figure  2, 
illustrates  the  major  features  of  the  system  which  allows  the 
scanning  of  spent  fuel  to  be  accomplished  in  the  MITR  spent 
fuel  storage  pit.   The  detector  dewar  is  positioned  on  the 
upper  movable  carriage  and  is  surrounded  by  lead  shielding  to 
reduce  background.   The  upper  movable  carriage  is  positioned 
on  the  lower  movable  carriage,  which  permits  scanning  in  two 
perpendicular  directions  in  the  horizontal  plane.   Extending 
from  the  upper  carriage  through  the  water  to  the  fuel  element 
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DETECTOR  DEWAR- 
LEAD  SHIELDING 


6  FEET  BELOW 
SURFACE  OF  WATER 


SPENT  FUEL7 
STORAGE  PIT 


GAMMA  RAY 
BEAM  TUB 5 


MITR  FUEL- 
ELEMENT 

POSITIONING  MECHANISM 


FIGURE  2    SCHEMATIC  DIAGRAM  of  APPARATUS  for  GAMMA  RAY 
SCANNING  of  SPENT  FUEL  ELEMENTS 
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is  a  1/2  in.  I.D.  air  filled  tube  which  acts  as  a  collimator 
and  allows  a  well-collimated  beam  of  gamma  rays  to  reach  the 
detector.   Additional  collimators  made  of  lead  with  diameters 
from  1/2  in.  to  1/8  in.  may  be  placed  in  the  beam  of  the  top 
carriage  to  reduce  the  count  rates  to  acceptable  levels  and 
provide  additional  shielding  against  background.   In  the  pre- 
sent experiment,  lead  sheets  were  also  placed  in  the  beam  to 
absorb  many  of  the  low  energy  gamma  rays  in  the  beam  while 
allowing  the  higher  energy  gamma  rays  to  pass  through  rela- 
tively unattenuated.   Count  times  of  80  to  160  minutes  were 
used  depending  on  the  count  rate. 

A  schematic  diagram  of  the  electronic  components  used 
together  with  important  settings  are  shown  in  Figure  3.   A 
list  of  these  components  is  shown  in  Table  10.   The  settings 
shown  in  Figure  3  resulted  in  an  energy  decrement  of  approxi- 
mately .66  keV/channel. 

3.3   Experimental  Results 

Identification  of  the  photopeaks  noted  in  the  spectra 
obtained  from  the  spent  fuel  elements  was  made  utilizing  the 
energy  calibration  of  Reference  (5) .   The  improved  resolu- 
tion of  the  Ge(li)  detectors  now  manufactured  at  M.I.T.  allow- 
ed better  resolution  of  the  spectra,  while  their  increased 
size  and  hence  greater  efficiency  made  shielding  against  back- 
ground more  difficult.   Orginally,  a  detector  with  an  active 
volume  of  30  cm  was  used,  but  even  with  4  inches  of  lead 
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Ge(Li)  DETECTOR 

PREAMPLIFIER 

BIAS  SUPPLY 
BIAS  =500  VOLTS 

AMPLIFIER 
COARSE  GAIN  =  8 
FINE  GAIN  =8.5 


4096  CHANNEL  ANALYZER 
DISCRIMINATOR  =1.7 


MEMORY 


READOUT  CONTROL 


SPECTRUM  STABILIZER 


OSCILLOSCOPE 


FIGURE  3    SCHEMATIC  DIAGRAM  of  ELECTRONICS  USED  for 
OBTAINING  GAMMA  RAY  SPECTRA  FROM  SPENT 
MITR  FUEL 
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TABLE  10:   ELECTRONIC  EQUIPMENT  DATA 


EQUIPMENT 


MANUFACTURER 


MODEL  NO 


Preamplifier 

Bias  Supply 

Amplifier 

4096  Channel  Analyzer 

Spectrum  Stabilizer 

Memory  Unit 

Readout  Control 


Canberra  Instruments 
Canberra  Instruments 
Sturrup 
Nuclear  Data 
Nuclear  Data 
Nuclear  Data 
Nuclear  Data 


1408C 

3001 

1417 

ND-161F 

ND-502 

ND-160M 

ND-160R 
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shielding  on  the  bottom  of  the  detector  and  2  inches  of  lead 
shielding  on  the  side  of  the  detector,  dead  times  of  30-35% 
were  noted  even  with  the  collimator  opening  also  shielded  by 
4  inches  of  lead.   Additional  shielding  could  not  be  added  as 
the  amount  of  shielding  material  was^ ^.limited  by  the  mech- 
anical strength  of  the  scanning  equipment.   This  larger  detector 

3 

was  replaced  with  a  smaller  detector  having  an  8   cm   volume. 

The  resulting  dead  times  were  found  to  be  on  the  order  of 
30%  with  the  collimator  tube  open  and  an  additional  1/8  in. 
aperture  collimator  inserted  immediately  below  the  detector 
when  element  7M9  was  placed  in  the  holder. 

3.3.1)   Fission  Product  Gamma  Ray  Spectra 

3.3.1a)   Fuel  Element  7M9 

The  spectrum  obtained  from  MITR  fuel  element  7M9 
after  a  cooling  period  of  79  days  is  shown  in  Figure  4.   The 
spectrum  was  obtained  with  a  1/8  in.  diameter  aperture  colli- 
mator for  a  counting  time  of  80  minutes.   In  addition  three 
1/4  in.  lead  sheets  were  placed  in  series  in  the  beam  to  pre- 
ferentially reduce  low  energy  gamma  rays  in  the  beam.   The 
spectrum  shows  resolved  peaks  for  gamma  rays  from  the  decay 
of  Nb-95,  Pr-144,  Cs-134  and  Cs-137.   Except  for  the  peaks  due 
to  Pr-144,  all  resolved  peaks  lie  between  500  and  800  keV,  and 
an  enlargement  of  this  portion  of  the  spectrum  is  shown  in 
Figure  5. 
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FIGURE  5    GAMMA  RAY  SPECTRUM  (600~800KEV)  OF  MITR 
FUEL  ELEMENT  7M9  AFTER  79  DAYS  COOLING 
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3. 3. lb)   Fuel  Element  5M2  7 

The  spectrum  from  MITR  fuel  element  5M27  after  a 
cooling  period  of  272  days  is  shown  in  Figure  6.   The  spec- 
trum was  obtained  using  a  1/8  in.  diameter  aperture   colli- 
mator for  a  counting  time  of  160  minutes.   One  1/4  in.  lead 
sheet  was  placed  in  the  beam  to  preferentially  reduce  the  low 
energy  gamma  rays  in  the  beam.   The  spectrum  shows  resolved 
peaks  for  gamma  rays  due  to  the  decay  of  Nb-95,  Zr-95,  Pr-144, 
Cs-134  and  Cs-137.   The  spectrum  between  600  key  and  800  keV 
is  shown  in  Figure  7. 

3.4   Calculations  of  Gamma  Ray  Intensities 

Intensities  of  gamma  rays  resolved  in  the  spectra  of 
Figures  4  and  6  were  calculated  by  finding  the  area,  express- 
ed in  terms  of  counts,  under  each  peak.   These  calculations 
were  performed  using  GAMANL,  Reference  (16) ,  a  computer  code 
which  identifies  peaks  in  spectra  and  determines  their  cen- 
ters and  areas.   GAMANL  utilizes  Fourier  transforms  to  smooth 
the  data,  which  eliminates  much  of  the  random  fluctuations  in 
the  spectra  without  affecting  the  resolution.   The  background 
is  then  calculated  and  subtracted  from  the  smoothed  data  and 
the  resulting  peaks  are  then  analyzed  for  peak  centers,  areas 
and  FWHM's,  among  other  things.   The  areas  are  calculated  in 
two  ways,  in  the  first  method  all  the  counts  are  added  be- 
tween successive  zeroes  in  the  smoothed  data  with  background 
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subtracted.   The  second  method  uses  a  Gaussian  fit  to  the 
peaks  and  calculates  the  area  using  the  height  of  the  peak  and 
its  standard  deviation. 

Results  for  the  spectra  shown  in  Figures  4  and  6  are 
shown  in  Tables  11  and  12  respectively. 

3.5   Conclusions 

Based  upon  the  results  obtained  in  the  work  reported  in 
this  chapter,  the  following  conclusions  can  be  drawn  pertinent 
to  the  projected  future  experiments  on  power  reactor  fuel: 

1)   The  improved  resolution  of  the  Ge(Li)  detectors  now 
manufactured  at  M.I.T.  and  the  4096  channel  and  analyzer  are 
evident  in  partial  resolution  of  the  Zr-95   Nb-95  doublet  at 
758  and  766  keV.   This  consideration  is  not  too  important  for 
cooling  times  on  the  order  of  one  year  or  less  since  the  peaks 
necessary  for  analyses  are  all  strong  peaks  and  tend  to  mask 
weaker  peaks  in  their  Compton  background.   After  several 
years  cooling,  however,  the  peaks  due  to  gamma  rays  from  Pr-144 
and  Cs-134  will  have  decreased  considerably  due  to  their  rela- 
tively short  half  lives  of  285  days  and  one  year  respectively. 
The  peaks  due  to  Nb-95  and  Zr-95  will  decrease  even  more 
rapidly  due  to  their  effective  half  life  of  65  days,  and  will 
be  a  relatively  small  peak  after  a  cooling  time  on  the  order 
of  1  1/2  years.   The  activity  of  fission  products  with  longer 
half. lives,  although  having  lower  fission  yields,  will  begin 
to  approach  the  activity  of  Zr-95,  Nb-95,  Pr-144  and  Cs-134, 
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and  the  resolution  of  the  system  becomes  considerably  more  im- 
portant. 

2)  In  the  future,  the  discriminator  setting  should  be 
set  higher  (e.g.  at  400  keV)  than  it  was  for  these  runs  (52 
keV)  to  help  reduce  counting  system  dead  times.   Because  it 
was  desirable  to  minimize  the  need  for  movement  of  the  detect- 
or, the  first  fuel  element  was  placed  in  the  holder  prior  to 
placing  the  dewar  containing  the  detector  on  the  top  carriage. 
This  made  it  impossible  to  make  background  runs  prior  to  ob- 
taining the  original  spectra,  and  peaks  noted  in  the  channels 
below  900  were  at  one  time  thought  to  be  significant.   Subse- 
quent background  runs  proved  this  original  assumption  to  be 
incorrect.   With  MITR  fuel  element  7M9  in  the  holder  but  with 
the  collimator  tube  removed  dead  times  on  the  order  of  15% 
were  noted  for  a  discriminator  setting  of  1.7,  corresponding 
to  a  cut  off  energy  of  52  keV.   For  a  discriminator  setting 

of  4.2,  corresponding  to  a  cut  off  energy  of  350  keV,  the  dead 
times  were  reduced  to  3%. 

3)  With  the  reduced  background  that  will  result  from 

the  less  active  spent  power  fuel  pins,  it  would  be  possible  to 

use  a  larger  detector,  which  with  its  improved  intrinsic 

efficiency  will  increase  the  number  of  counts  measured, 

especially  for  the  higher  energy  gamma  rays.   The  use  of  a 

larger  detector  should  also  help  reduce  the  Compton  edges 

somewhat  by  enabling  more  gamma  rays  to  deposit  their  full 

energy. 
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4)   Calibration  of  the  system  for  the  double  escape  peak 
of  the  2186  kev  Pr-114  gamma  ray  should  also  improve  the  deter- 
mination of  the  Pr-144  content  of  the  spent  fuel.   However,  it 
should  be  noted  the  use  of  a  larger  detector  may  decrease  the 
probability  of  both  gamma  rays  escaping  from  the  detector  after 
a  Pair  production  event  and  thus  decrease  the  size  of  the 
double  escape  peak. 

In  short,  it  is  concluded  that  the  proposed  power  fuel 
experiments  using  the  apparatus  described  in  this  chapter  are 
quite  feasible.   Sufficient  signal  strength  is  easily  obtain- 
ed, only  l/16th  of  the  full  collimator  bore  was  used  in  the 
present  scoping  studies.   Room  background  can  be  reduced  to 
acceptable  levels  by  proper  choice  of  detector  size,  lead 
filter  thickness  and  shielding.   Dead  time  can  be  reduced  by 

Using  a  sufficiently  high  low  level  discriminator  cut  off. 
The  key  limitation  is  really  in  the  Compton  background  pro- 
duced in  the  detector. 

So  far,  then,  we  have  investigated  and  established  the 
potential  ability  to  handle  spent  fuel  and  carry  out  a  spec- 
trometer scan  with  available  apparatus.   In  the  next  chapter 
a  projection  will  be  made  concerning  the  suitability  of  the 
data  obtainable  with  this  apparatus  for  the  fuel  in  question. 
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IV:   CALCULATION  OF  GAMMA  RAY  INTENSITIES  EXPECTED  IN 
SPENT  FUEL. 

4 . 1  Introduction 

In  order  to  determine  the  feasibility  of  studying  the 
spent  fuel  made  available  by  the  AEC  it  was  necessary  to 
determine  the  expected  intensities  of  those  gamma  rays  use- 
ful for  analysis.   It  was  also  necessary  to  determine  whether 
they  would  be  resolvable  for  a  cooling  time  of  approximately 
2.75  years.   In  addition  to  calculating  the  gamma  rays  emitted 
by  fission  products  of  U-235  fissions,  an  estimate  was  made  of 
the  expected  gamma  rays  from  fission  products  of  Pu-239. 

4 . 2  Calculation  of  Expected  U-235  Fission  Product  Gamma 
Ray  Intensities. 

Calculation  of  the  expected  intensities  for  gamma  rays 
from  U-235  fission  products  was  performed  by  utilizing  the 
calculated  peak  area  (expressed  as  counts)  obtained  in  Chapter 
3.   The  data  obtained  for  MITR  fuel  element  5M27  was  extra- 
polated, using  the  familiar  experimental  decay  relationships, 
to  a  2.75  year  cooling  period,  which  is  the  expected  cooling 
period  for  the  Dresden  fuel  pins.   The  calculated  intensities 
are  shown  in  Table  13. 

In  order  to  check  the  validity  of  extrapolating  the 
intensity  of  a  gamma  ray,  as  measured  by  the  area  under  its 
peak,  on  the  gamma  ray  spectra,  a  parallel  calculation  was 
performed  extrapolating  the  measured  intensity  for  the  gamma 
rays  of  7M9  (cooling  time  =  2.75  months)  to  a  cooling  period  of  9 
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TABLE  13 


INTENSITY  OF  U-2  35  FISSION  PRODUCT  GAMMA  RAYS 


AFTER  2.75  YEARS  COOLING 


FISSION 
PRODUCT 


ENERGY 
(keV) 


PEAK 
AREA 


Pr-144 


2186 
1490 
1164 
696 


322.8 
174.5 
992.4 
1,102. 


Cs-134 

796 

4,129.3 

605 

4,327.4 

Nb-95 

766 

93.3 

Zr-95 

758 

32.9 

724 

19.9 

Cs-137 

662 

29,565.7 

Rh-106 

623 

544.0 

511 

1280.0 
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months.   Corrections  were  made  for  the  differences  in  counting 
time,  thickness  of  lead  sheets  in  the  beam,  and  burnup.   These 
calculated  intensities  are  compared  with  the  measured  inten- 
sity (peak  area)  for  the  gamma  rays  of  5M2  7  (cooling  time  =  9 
months)  in  Table  14. 

The  calculated  intensities  compare  well  with  the  measur- 
ed intensities  with  the  exception  of  Cs-137.   No  reason  for 
the  discrepancy  could  be  discovered  since  both  5M2  7  and  7M9 
had  similar  irradiation  histories  and  were  used  in  similar 
positions  in  the  MITR  core. 

In  order  to  determine  if  the  values  of  the  calculated 
intensities  listed  in  Table  13  were  sufficiently  strong  to  be 
resolved,  their  relative  intensities  were  analyzed.   This 
method  of  analysis  was  chosen  because  with  spent  fuel  gamma 
spectra  the  difficulty  is  not  one  of  obtaining  sufficient 
count  rate  but  rather  of  resolving  lower  energy  gamma  rays  in 
the  presence  of  Compton  background  from  higher  energy  gamma 
rays.   Since  the  height  of  the  Compton  edge  and  continuum  are 
proportioned  to  the  size  of  the  peak  it  was  felt  that  the 
ratio  of  a  gamma  ray's  peak  area  to  that  of  the  sum  of  the 
peak  areas  for  higher  energy  gamma  rays  would  yield  sufficient 
information  to  predict  which  gamma  rays  would  be  resolvable. 
The  results  of  this  analysis  are  shown  in  Table  15.   For  com- 
parison the  results  of  the  GAMANL  (16)  calculation  for  fuel 
element  5M2  7  were  analyzed  in  this  manner  and  the  results  are 
listed  in  Column  3  of  Table  15.   From  this  it  was  concluded 
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TABLE  14:   CALCULATED  AND  MEASURED  INTENSITY  OF  U-2  35  FISSION 
PRODUCT  GAMMA  RAYS  AFTER  9  MONTHS  COOLING. 


ENERGY 

CALCULATED 

MEASURED 

(keV) 

PEAK. 

PEAK 

AREA 

AREA 

2186 

1,836.9 

1914.0 

1490 

1,105.3 

1034.5 

1164 

6,904.4 

5884.4 

796 

8,880.6 

8003.9 

766 

178,525.2 

228,593.6 

758 

69,640.8 

80,577.3 

724 

46,819.0 

48,988.5 

696 

6,980.4 

6,535.3 

662 

53,573.6 

30,966.5 
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TABLE  15 


CALCULATED  VALUES  OF  RELATIVE  INTENSITY  FOR  U-2  35 


FISSION  PRODUCTS  AFTER  2.75  YEARS  COOLING. 


ENERGY 

CALCULATED 

MEASURED 

(keV) 

RELATIVE 

RELATIVE 

INTENSITY 

INTENSITY 

(9  mos.  cooling) 

2186 

_ 

_ 

1490 

.54 

.54 

1164 

1,99 

1.99 

796 

8.30 

2.72 

766 

.02 

20.70 

758 

.0069 

.  337 

724 

.0042 

.153 

696 

.227 

.017 

662 

5.03 

.082 

623 

.015 

.008 

605 

.106 

.020 

511 

.030 

.017 
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that  all  gamma  rays  were  capable  of  being  resolved  by  use  of 
GAMANL  (16).   However,  it  is  felt  that  the  low  absolute  in- 
tensity of  the  peaks  at  766,  758  and  724  keV  would  not  lead  to 
reproducible  results  in  the  calculation  of  the  areas  under 
the  peaks. 

4 . 3   Calculation  of  Expected  Pu-239  Fission  Product 

4.3.1   Gamma  Ray  Intensities 

In  order  to  obtain  an  estimate  of  the  number  of 
Pu-2  39  fissions  which  have  occurred  in  the  fuel  pins  an  energy 
balance  was  performed  using  power  history  data  supplied  by 
General  Electric  for  the  fuel  in  question.   It  was  assumed 
that  the  energy  was  released  either  by  thermal  fission  of  U- 
235  or  Pu-239,  i.e. : 

E (TOTAL)  =  E (Pu-239)  +  E(U-235) 

Assuming  200  Mev  are  released  per  fission: 

2  3 
,  ,.  ,.     ,„    6.023  x  10    atoms     gm  mole        1  fission 

gm  -  mole         2  35  gm  U-2  35      1  atom 


200  Mev    577h       fission  U-235           watt  -  sec 
x  x  u     x  

fission    683b  fission  +  capture  U-235    6.25  x  1012  Mev 

x  day       x     MW     =  15.3(6.023)  (577)  (200)  =    12    2   MWD 

8.64  x  104  sec    106  watt  =  235  (683)  (8. 64)  (6 . 25) 


Thus: 

E (Pu-239)  =  E (TOTAL)  -  E  (U-2  35) 
-50- 


=  17.08  -  12.20 

=  5 . 0  8  M WD 

m,     ,  .    ,.   fissions  in  Pu-239    5.08    .-, 
The  ratio  of   =  — - =.416 

fissions  in  U-235     12.2 

Since  approximately  42%  of  the  fissions  are  in  Pu-239 
there  is  a  possibility  of  resolving  Pu-239  fission  product 
gamma  rays  in  the  spectrum.   The  value  of  42%  is  in  good 
agreement  with  the  values  predicted  using  cross  section  data 
as  reported  in  Appendix  A. 

From  Table  8  it  is  seen  that  the  fission  yield  of  Rh-106 
-  Ru-106   is  increased  from  0.4%  for  thermal  fission  of  U-235 
to  4.4%  for  thermal  fission  of  Pu-2  39.   Table  13  shows  the 
predicted  peak  area  for  the  623  keV  peak  of   Rh-106  is  544.0 
counts  for  fission  of  U-235.   The  peak  area  for  fission  of 
Pu-239  would  be 

iii*  x  544.0  =  5440.0  counts 
.4% 

The  superimposed  peak  area  for  the  fuel  pins  for  the  623  keV 

gamma  ray  would  then  be 

.416(5440.0)  +  .584(544.0)  =  2574.0  counts 

This  is  an  appreciable  percentage  of  the  peak  area  for  Cs-137 
and  thus  should  be  resolvable. 

4.4   Conclusions 

A  sufficient  number  of  gamma  ray  peaks  should  be  resolv- 
able and  statistically  significant  in  the  spectra  to  enable 
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nondestructive  testing  to  be  performed.   Only  one  gamma  ray 
need  be  resolved  for  each  fission  product.   Even  though  the 
gamma  rays  of  Nb-95  and  Zr-95  will  not  be  resolvable,  the 
gamma  rays  of  Rh-106  at  62  3  keV  should  be  resolvable  and  thus 
the  method  of  analysis  utilizing  three  fission  product  activ- 
ity ratios  recommended  by  Sovka  (5)  should  be  possible. 
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V:   CONCLUSIONS 

The  major  conclusions  reached  as  a  result  of  the  work  re- 
ported in  this  thesis  are  as  follows: 

1)  The  calculated  dose  rate  for  the  Dresden  power  re- 
actor fuel  pins  is  well  below  that  presently  experi- 
enced from  spent  MITR  fuel  elements,  and  thus  pre- 
sents no  special  handling  problems. 

2)  Sufficient  spectral  information  should  be  obtained 
from  the  spent  Dresden  fuel  pins  to  enable  non-destruc- 
tive analysis  to  be  performed.   In  fact,  the  increased 
energy  resolution  offered  by  present  day  Ge  (Li)  detect- 
ors and  the  40  9  6  channel  analyzer  coupled  with  the  im- 
proved analytic  capability  of  the  GAMANL  code  (16) 
allow  improved  determination  of  the  actual  gamma  ray 
intensities  from  spent  fuel  spectra  as  compared  to 
previous  studies  on  MITR  fuel  by  Sovka  (5) . 

3)  Although  Nb-95  and  Zr-95  gamma  rays  will  probably  not 
be  resolved  in  the  subject  Dresden  fuel,  the  expected 
relative  increase  in  Ru-106  and  Rh-106  will  permit  sub- 
stitution of  this  latter  pair  in  the  analytic  procedure. 

4)  A  method  of  Compton  suppression  would  be  extremely 
valuable  because  of  the  intense  Compton  background  pre- 
sent between  600  and  800  keV,  the  energy  range  for  all 
but  one  of  the  gamma  rays  required  for  the  method  of 
analysis  presented  by  Sovka  (5).   In  addition,  Compton 
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4)   cont. 

suppression  might  enable  weaker  intensity  gamma 
rays  potentially  useful  for  analysis  to  be  re- 
solved.  The  methods  utilizing  anti-coincidence 
would  be  extremely  difficult  to  apply  because 
of  volume  and  weight  limitations  imposed  by  the 
scanning  equipment.   A  duode  type  Ge(Li)  detect- 
or (17,18),  however,  would  be  ideally  suited  to 
scanning  of  spent  fuel  since  it  requires  no 
additional  space  or  volume  than  that  needed  for 
a  standard  Ge(Li)  detector. 
As  a  result  of  this  work,  it  can  be  said  that  all  necessary 
groundwork  for  spent  power  reactor  fuel  pin  analysis  has  been 
completed  to  the  point  where  shipping  and  receipt  of  the  fuel, 
and  initiation  of  scanning  experiments  is  now  warranted. 
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APPENDIX  A 

A  parallel  calculation  of  the  ratio  of  the  number  of  fission 
in  Pu-2  39  to  the  number  of  fissions  in  U-2  35  was  carried  out  us- 
ing one  group  theory  utilizing  effective  neutron  cross  sections 
to  account  for  the  fission  caused  by  epithermal  neutrons  and  to 
account  for  all  neutron  captures  in  U-238. 

In  order  to  obtain  an  expression  for  the  number  of  Pu-239 
fissions  we  must  first  obtain  an  expression  for  the  number  of 
U-238  atoms  converted  to  Pu-239  atoms  as  a  function  of  time. 
For  the  one  group  model  we  can  write: 


dN 

at 


28 


-(T  286u2i 


(1) 


where : 


N 


<£,Y 


28 


28 


=  no.  of  U-2  38  atoms  at  time  t 

=  effective  microscopic  capture  cross  section 
for  U-238 

=  one  group  neutron  flux 


Solving  Eq.  (1) 


N28  =  N28 
o 


exp("°rTY   ^fc) 


(2) 


where : 


2  8 
N    =  initial  no.  of  U-238  atoms 
o 


For  Pu-2  39  atoms  we  can  write 


§f  =  ^,Y2>28  -  ^  *49 


(3) 
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where  : 


49 
N        =   no.    of   Pu-2  39    atoms    at    time    t 

49 
(p         =   effective   microscopic    cross    section    for   Pu- 
239 


2  8 
Substituting   Eq.     (2)     for  N         in   Eq.     (3) 


49 

dN  ^_     28  /  ^,28 

dt     =  <5TfY    9  No 


exp (- 0 


28 


^  t) 


49       49 
-  (j-         N4y  (4) 

d 


Solving   Eq.     (4)    yields 


49  49 

N        =   N 
o 


exp(-(7^4Y  t)    +  N    1       <L*il 


28  \ 


28  49 

^n,Y        +  ^a         /   _ 


^P^^y28^^" 


exp(-(J^      0t) 


49  _49  28  -2.3 

Since   N  =0.    <j  ~7P' CT~  and  exp(-    (j  t)  sss  1    for  normal 

fuel   burn up    the    above   expression    can   be    approximated  by 


N         =   N 
o 


f  28\ 

(Tn",y 


*~a 


49 


,—-49  / 
l-exp(-  <Ja<pt) 


(5) 


For  the  time  ratio  of  Pu-239  fissions  we  can  write 


dF 


49 


-49  J  .,49 


K  =  «T*" 


(6) 


49 
Substituting  Eq.  (5)  for  N    in  Eq. (6)  we  have 


dF 

at 


49 


49i  ..  2 


-  tf  Y  % 


28 


tf-"49  / 
a 
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49 


l-exp(-  (T  (f>   t) 


Integration  yields 


49  -49jxt21 

F         =     q       0NQ 


£■, 


28 


Y 


(T 
a 


49 


t- 


l-exp((f49^  t) 


49 


<V 


(7) 


Expanding    the   exponential   term   of    Eq.     (7)    in    a   Taylor   Series 
yields 


.49 


0No 


/  28 

(fn,Y 


1 


<t 


49 


49  /  49/2 

t-    1-1    +    (T      0  t-l/2((T^  t)     +... 


\ 


,—49/ 


/ 


Ignoring   all    terms    above    the    t      term  yields 


/  281 


49  _.49j-„t2  8 


Gn, 


Y 


49,2 


i<S 


49 


(1/2    Gl     (ft  t*) 


(8) 


For  the  number  of  U-2  35  atoms  we  may  write 


dN 
dt 


25 


-  e:254  n25 


(9) 


where 


.25 


N    =  no.  of  U-235  atoms  at  time  t 
.25 

U-235 


Q~Z~      =  effective  microscopic  cross  section  for 


Solving  Eq.  (9) 


N25  =  N^5  exp(-  r25Cp  t) 


(10) 
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The  rate  of  U-235  fissions  is  given  by 


dF 


25 


25  /  25 


-r"4> 


N 


(ID 


where 


.25 


F   -  total  no.  of  U-235  fissions  from  time  t  =0 
Substituting  Eq.  (10)  into  Eq.  (11) 


dF 


25 


25 


=  fjl   VNo     (l-exp(-(T25^t))/(r 


.2  5 


(12) 


The    ratio   of   Pu-239    fissions    to   U-235    fissions    using  Eq.     (8) 
and  Eq.     (12)    becomes 


F49  [.ff 

175"  " 


49\ 


L.    28  2  3\/M  /oV_49   /  ^2.  ,^-25y    . 


\ 


r-49 


"o"    tf  ' 


25         25 
1-     (N/D/No       ) 


/ 


67 


49 


N 


28 


**, 


28\ 


T 


,25         VT25        ^-25 
\   o  '  \       f        / 


^2V>2 


(13) 


In   order   to   eliminate    flux   and   time    dependence   of   Eq.     (13 
return   to   Eq.     (10)    and   take    logarithms    of   both   sides 


.     /M25  /1T25,  ,-2  5  /   . 

In  (N      /Nq    )    =    $2      (f)  t 


Solving    for  <fit  yields 

0t   = 


25  25       25 

(i/(T^    )  murvN*5; 
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Substituting  for   t  in  Eq.  (8)  yields 


.49 


.25 


1/2 


67 


49 


^25, 


/  28\/ 


SI 


w 


N 


28 
o 


25M  25 


N 


25 


In 


'n25 

o 


,T25 
N    I 


^2 


(14) 


Rearranging. yields 


.49 


I  rr*9\L28u    28\ 


=    1/2 


F 


25 


£ 


(T      N 


ni 


^25T7  25 
0    N 
a   o  / 


H»2X"  / 


/.  25\ 


In 


N 
o 


N 


25 


O  Q     OQ        O^     OCT 

The  ratio  C~     N    /  0~  N     is  defined  as  the  initial  conver- 
a   o  /   a   o 

sion  ratio  (ICR)  and  is  equal  to  the  initial  number  of  neutron 
absorptions  in  U-238  divided  by  the  number  of  neutron  absorp- 
tions in  U-2  35.   Thus  the  final  working  equation  becomes 

2 


,49 
"25" 


=  1/2 


<ff 


49 


ICR 


<rf25J   (i-n25/n025> 


In 


N 


25 


o 


N 


25 


(15) 


In  order  to  apply  Eq.  (15)  we  need  to  specify  the  effec- 

49      ,—49 
tive  cross  sections  (JI   and  <j   .   To  account  for  both  ther- 

f        a 

mal  and  epithermal  events  the  following  approximate  formula- 
tion can  be  employed: 


cr   -  ck  + 


^25  WRI 


25 


RI 


J  f 


where:   0~~  =  effective  microscopic  cross  section 
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(TTV  =  thermal  microscopic  cross  section 
(jT(.  =  measured  ratio  of  epithermal  to  thermal  fissions 
RI  =  infinite  dilution  resonance  integral 


Using  tabulated  values  for  (T7h(19)  and  RI(20)  and  a  value  of 

O^r    =  0.144,  a  typical  value  for  the  lattice  in  question,  we 

49  _25 

obtain   (TT    =  826  barns  and  cr;   =  660  barns.   Using  these 
t  * 

49         25 
values  for   (T.    and  G7   ,  a  typical  value  for  ICR  of  0.5, 

25       25 

and  the  values  of  N     and  N    from  Table  1  we  obtain  from 

o 

Eq.  (14)  a  value  of  F49/F25  =  .392. 
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